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1 Introduction 



During the past decades, there has been an increasing interest in the study of classical and quantum 
integrability of Heisenberg ferromagnet (HM) model [T]-[15j. The Heisenberg ferromagnet (HM) 
model based on Hermitian symmetric spaces has been studied in [Il]-[I1]. The integrability of 
the HM model based on SU{2) via inverse scattering method is presented in [2]- [3] and its SU{n) 
generalization is studied in [4J . The integrability of the GHM model based on the general linear Lie 
group GL{n) via Lax formalism has been investigated in [1]. In this paper we present the Darboux 
transformation of the GHM model based on general linear group GL{n) with Lie algebra gl(n) 
and calculate multi-soliton solutions in term of quasideterminants. We also establish the relation 
between the Darboux transformation and the well-known dressing method [16]. In the last section, 
we discuss the model based SU (n) and calculate an explicit expression of the single-soliton solution 
of the HM model based on the Lie group SU{2) using Darboux transformation. 

The Hamiltonian of the GHM model is defined by [T] 



with "T" is transpose and U{x,t) is a matrix- valued function which takes values in the Lie algebra 
gl(n) of the general linear group GL{n). The corresponding equation of motion can be expressed 
as 




(LI) 



dtU = {H,d^u]. 



(L2) 



The above equation (jl.2p can be written as 



dtU=[u,dlu] 



(1-3) 



where dx = qz and dt = ^. Let us assume that U{x,t) is diagonizable, i.e.. 



U = gTg 



1 



(L4) 



where g G GL{n) is matrix function of (x, t) and T is a n x n constant matrix 

/ci ••• ••• o\ 
ci ■■■ 0---00 







ci 







(L5) 



T = 







C2 











C2 







\ ••• ••• 0c2/ 
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where 1 < p < n and ci , C2 € K (or C) . From equations (jl.4p and (jl.Sp , we have 

[U,[U,[U,x]]]=c'' [U,x], 
for an arbitrary matrix function x and c = ci — C2 7^ 0. Since 

imphes 

[U,[U,U^]]=c^U,, 

The equation of motion (jl.3p can also be written as the zero-curvature condition i.e., 



(1.6) 



(1.7) 



1 



-u,dt 



1 



[u, u^] 



0. 



:ii 



(l-A)"'"^ (1_A)2 (i-A) 
The above zero-curvature condition (|1.9p is equivalent to the compatibility condition of the following 
Lax pair 

1 



(1-A) 



U{x,t)'^{x,t;X) 



1 



[/,[/,] ^{x,t;X) 



(1.10) 
(1.11) 



.(1-A)2 (1-A) 

where A is a real (or complex) parameter and ^ is an invertible n x n matrix-valued function 
belonging to GL{n). 

In the next section, we define the Darboux transformation on matrix solutions ^ of the Lax pair 
()1.10p - ()l.lip . To write down the explicit expressions for matrix solutions of the GHM model, we 
will use the notion of quasideterminant introduced by Gelfand and Retakh |17j-|21j. 



Let X be an n X n matrix over a ring R (noncommutative, in general). For any 1 < i, j < n, 
let Tj be the ith row and Cj be the jth column of X. There exist n? quasideterminants denoted by 
\X\ij for i, j = 1, . . . , n and are defined by 



\X 



(1.12) 



where Xij is the ijth entry of X, represents the ith row of X without the jih. entry, represents 
the jih. column of X without the iih entry and X^^ is the submatrix of X obtained by removing 
from X the ith row and the jih. column. The quasideterminats are also denoted by the following 
notation. If the ring R is commutative i.e. the entries of the matrix X all commute, then 



\X\ 



: deiX 
deiX^j 



(1.13) 
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For a detailed account of quasideterminants and their properties see e.g. [I7|-[2I]. In this paper, we 
will consider only quasideterminants that are expanded about an nx n matrix over a commutative 
ring. Let 

(c «)• 

be a block decomposition of any K x K matrix where the matrix D is n x n and A is invertible. 
The ring R in this case is the (noncommutative) ring of n x n matrices over another commutative 
ring. The quasideterminant oi K x K matrix expanded about the n x n matrix D is defined by 



A 
C 



B 



D 



D-CA-^B 



(1.14) 



The quasideterminants have found various applications in the theory of integrable systems, where 
the multisoliton solutions of various noncommutative integrable systems are expressed in terms of 
quisideterminants (see e.g. [22]- [30]). 



2 Darboux transformation 



The Darboux transformation is one of the well-known method of obtaining multi-soliton solutions 
of many integrable models [3T]-[33|- We define the Darboux transformation on the matrix solutions 
of the Lax pair (jl.lOp - ijl.lip . in terms of an n x n matrix Z)(x, t. A), called the Darboux matrix. For 
a general discussion on Darboux matrix approach see e.g. [34]- [39]. The Darboux matrix relates 
the two matrix solutions of the Lax pair (jl.lOp - (jl.lip . in such a way that the Lax pair is covariant 
under the Darboux transformation. The one-fold Darboux transformation on the matrix solution 
of the Lax pair (frTO]) - (frTT]) is defined by 

^ [1] (x, t; A) = D{x, t, A)^'(x, t; A), (2.1) 

where D(x,t,X) is the Darboux matrix. For our case, we can make the following ansatz 

D{x,t,X) = XI - M{x,t), (2.2) 

where M{x, f) is an n x n matrix function and / is an n x n identity matrix. The new solution 
^' [1] (x, t] A) satisfies the following Lax pair, i.e. 

d,^[l\{x,t-X) = ^C/[l]iI/[l](x,t;A), (2.3) 
dt^[l]{x,t-X) = (^^^-^C/[l] + ^[C/[l],C/,[l]]^xI/[l](x,i;A), (2.4) 
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where U [1] satisfies the equation of motion (jl.Sp . By operating dx and dt on equation (j2.ip and 
equating the coefficients of different powers of A, we get the following transformation on the matrix 
field U 

U[l] = U + Mx, (2.5) 

and the following conditions which M is required to satisfy 

Mx{I-M) = [C/,M], (2.6) 
Mt {I - Mf = [c^U + [U, C/^] , M] + M [U, Ux] M - [U, U^] M\ (2.7) 

One can solve equations (I2.6|) - p.7p to obtain an explicit expression for the matrix function M{x, t). 
An explicit expression for M{x, t) can be found as follows. 

Let us take n distinct real (or complex) constant parameters Ai,--- ,A„(7^ 1). Also take n 
constant column vectors €1,62,- ■■ , e„ and construct an invertible non-degenerate nx n matrix 
function Q(x,t) 

@{x, t) = (^(Ai)ei, • • • , ^(A„)e„) = {6^, ■ ■ ■ ,9^) . (2.8) 

Each column 6i = ^'(Aj)ej in the matrix is a column solution of the Lax pair ()1.10p - ()l.ll|) when 
A = Aj and i = 1, 2, . . . , n i.e. 

dx9i = r^^^*' (2-9) 

Let us take an n x n invertible diagonal matrix with entries being eigenvalues Aj corresponding to 
the eigenvectors 6i 

A = diag(Ai,...,A„). (2.11) 
The n X n matrix generalization of the Lax pair ()2.9p - ()2.10p will be 

d^e = c/e(/-A)"\ (2.12) 
dto, = c^ue{i-Ar^ + [u,u.,]@{i-AyK (2.13) 

The nxn matrix is a particular matrix solution of the Lax pair (I2.9p - (12.10p with A being a matrix 
of particular eigenvalues. In terms of particular matrix solution of the Lax pair (j2.9p - (|2.10p . we 
make the following choice of the matrix M{x, t) 

M{x,t) = QAQ^'^. (2.14) 



5 



Our next step is to check that equation (j2.14p is a solution of equations (j2.6p - ()2.7p . In order to 
show this, we first operate dx on equation (I2.14p to get 

dxM = 5,(eAe-i), 

= (9,e)AG-i + GAa,(e'i), 

= UQ{I - K)-^KQ-^ - QKQ-^Ue{I - A)-^e^\ 

= -U + {I - M)U{I - Uy^ , (2.15) 
which is the equation (j2.6p . Similarly operate dt on (j2.14p . we get 

dtM = dt{eAe-^) 

= (dtQ) AQ-^ + eQAdt{@-^) 

= (^c^ue {I - A)-2 + [u, Ux] e (/ - A)-^) AQ-^ - 

QAQ-^ (c^UQ (I - A)-^ + [U, Ux] G (/ - A)"^) Q-\ (2.16) 



which is equation (j2.7p . This shows that the choice (j2.14p of the matrix M satisfies the equations 
(j2.6p - (j2.7p . In other words we can say that if the collection (^,U) is a solution of the Lax pair 
(frT0]) - (frTT1) and the matrix M is defined by (EH]), then (^'[1], C/[l]) defined by I^Ji) and (|23]) 
respectively, is also a solution of the same Lax pair. Therefore we say that 

^[1] = (A/-eAG"^)^, 

U[l] = {I- QAQ-^) U{I- OAQ-^) , 

is the required Darboux transformation on the solution ^ to the Lax pair p.lUp - p.lip and U to 
the equation of motion (jl.Sp respectively. 



3 Quasideterminant solutions 



We have shown that the matrix M = GAG ^ satisfies the conditions (j2.6p - (|2.7p . Therefore, the 
one-fold Darboux transformation (|2.ip can also be written in terms of quasideterments as 



^-[1] = L>(x,t;A)^' = (A/-GiAiGf^)^', 
Gi ^' 



GiAi 



(3.1) 



6 



The above equation defines the Darboux transformation on the matrix solution ^ of the Lax pair 
(|1.10P - (|l.lip . The corresponding one-fold Darboux transformation on the matrix field U is 



U[l] 



Gi / 61/ 

91 (/-Ai) 91 (/-Ai) [0] 



-1 



We write two-fold Darboux transformation on ^ as 



Z)(x,t;A)^[l] = A^[l] -92[l]A292Ml]^'[l] 
A (A/ - 9iAi9^^) ^ - 

(92A2 - 9iAi9r^92) A2 (92A2 - 9iAi9r^92)"^ (A/ - 9iAi9r^) ^, 

9i 92 
9iAi 92A2 A'f 

9iA2 92A2 



A"^ 



Similarly the expression for two-fold Darboux transformation on the matrix field U as 
U[2] = 92[l](/-A2)92^[l][/[l](92[l](/-A2)92^[l])-\ 



, -1 



(92A2 - 9iAi9ri92) (I - A2) (92A2 - 9iAi9f ^92 
(/ - 9iAi9^^) U{I- 9iAi9^^)~^ X 
92A2 - 9iAi9ri92) (/ - A2) (92A2 - 9iAi9r^92)" 



(3.2) 



(3.3) 



Gi 92 / 

9i (/ - Ai) 92 (/ - A2) 
9i(/-Ai)2 92(/-A2)' 



xU X 



9i 82 / 

9i(/-Ai) 92(I-A2) 
9i(/-Ai)2 92(I-A2)' [0] 



-1 



(3.4) 



The result can be generalized to obtain A^-fold Darboux transformation on matrix solution ^ as 



9i 


92 • 


• 9Ar 




9iAi 


92A2 • 


• 9^Ajv 




9iA? 


92A2 . 


• 9^A2, 


A^M/ 



9iAf 92A^ 



A^^ 



(3.5) 
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Similarly the expression for U[N] is 



U[N] 



ei(/-Ai: 



62 

62 (/ - A2) 



@i{i-Aif e2(/-A2)' 



ei(i-Ai)^ e2(/-A2) 

ei 62 

ei(/-Ai) e2(/-A2) 

Gi(/-Ai)2 62(1 -Aaf 



N 



ei(/-Ai)^ 62 (/-A2)^ 



On I 
On {I -An) 
e^v (/ - A,v)' 



On (/ - A^) 



N 







xU X 



Qn I 
On {I - An) 
Qn {I - A^)' 



e^(/-A;v)^ [0] 



(3.6) 



We now relate the quasideterminant solutions of GHM with the solutions obtained by dressing 
method and the inverse scattering method. For this purpose, we proceed as follows. From the 
definition of the matrix M, we have 



Me 



BA. 



(3.7) 



Let 6i and Oj be the column solutions of the Lax pair (jl.lOp - ()l.lip when A = Aj and A = Aj 
respectively i.e. 



MOi = XiOi, i = l,2,...,p 

MOj = Xj9j. j =p + l,p + 2,...,n 

Now we take Xi = fi and Xj = fi, we may write the matrix M as 



(3.5 



(3.9) 



where P is the hermitian projector i.e. P^ = P. The projector P satisfies P^ = P and P-^ = 1 — P. 
The projector P is hermitian projection on a complex space and P"*" as projection on orthogonal 
space. Now equation (j3.9p can also written as 



M 



{fi- H)P + fil, 



(3.10) 



where the hermitian projector can be expressed as 



(3.11) 
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The one-fold Darboux transformation (j3.ip on the matrix solution ^ can also be expressed in terms 
of projector P as 

^[1] = P(x,t;A)^ = ^— ^, (3.12) 

where I?(x,t;A) is the rescaled Darboux-dressing function i.e. 'D{x,t;X) = (A — /i) ^ D{x,t;\). 
Similarly the A^-fold Darboux transformation (|3.5p on the matrix solution ^ can also be written as 
(take P[l] = P) 

N~l . _ _ . 

^[N] = W [l-l^Ii^^J^P[N -k]]^. (3.13) 
fc=0 ^ ^ " ^"^-^ ^ 

Now we can express the A^-fold Darboux transformation p.6p on the matrix field U can be written 
as 

N-l . _ X N-l , _ . 

\J\N\ = TT (l- ^^~^-^^-^ pyN-k\\uW (/-^1^P[/]], (3.14) 

and hermitian projector is defined as 



m = e^\k\(e''MM^)' (3.15) 

The expressions (|3.13p and (j3.14p can also be written as sum of K terms |27] 

^[iV] = (3.16) 

fc=o ^ J 

and 
where 

Af-l , 
i=0 



4 The explicit solutions of the GHM model 



In this section we calculate explicit expression of soliton solution. First of all we will study GHM 
model based on SU{n). In this case the spin function U takes values in the Lie algebra su(n) so 
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that one can decompose the spin function into components U = U"'T"', and T", a = 1,2, ... ,n are 
anti-hermitian n x n matrices with normahzation Tr (T"T^) = ^5'^^ and are the generators of the 
SU{n) in the fundamental representation satisfying the algebra 

jabcj^c^ (4.1) 



rj-ia rpb 



where f"'^'^ are the structure constants of the Lie algebra su(n). For any X G su(n), we write 
X = X^T" and U" = -2Tr([/r'='). 

The matrix-field U belongs to the Lie algebra su(n) of the Lie group SU{n) therefore 

C/t = -U, Tt{U) = 0. (4.2) 



The equations (j2.ip - ()2.2p and (j2.5p define a Darboux transformation for the GHM model based on 
the Lie group SU{n). The new solution of the equation of motion p.3p U[l\ must be su(n) valued 
i.e. 

= -[/[!], Tr([/[1]) = 0, (4.3) 

therefore, we have the following conditions on the matrix M 

= -M, Tr(M) = 0. (4.4) 

In other words we want to make specific M to satisfy the (j4.4p . This can be achieved if we choose 
the particular solutions 0, at A = Aj, let us first calculate 

= {i-\iyU\u^e, + {i-\,)-^e\uej, (4.5) 

using equation (j4.2p the above equation (j4.5p becomes 

dje\e,) = 0, (4.6) 



when Aj ^ Xj (i.e. Aj = Xj). Similarly we can check 

dt(e\e,) = 0. (4.7) 



From the definition of the matrix M, we have 

e\{M^ + M)ej = {Xi + Xj)ej6j, (4., 



10 



when Aj 7^ Xj then the above expression (j4.8p implies 



9le, = 0. 



(4.9) 



The column vectors 9i are linearly independent and the equation (14. 9p holds everywhere. 
For the HM model based on SU{n), the constant matrix (jl.Sp becomes 



/2 









\ 








\ 




2 

n 



2 

n 











(4.10) 



Then C/^ becomes 



(4.11) 

These are the constraints given in ref. [4J . For the construction of explicit soliton solution for the 
SU{n) HM model, we construct the matrix M by defining a Hermitian projector P. For this case, 
we take the seed solution to be 



Un = U = -i 



( 



V 



ai 



(4.12) 



where ai are real constants and 'Y^=\ — 0- The corresponding solution of the Lax pair is expressed 
in block diagonal matrix 

' W,,{\) O \ 

O Wn-y{\) ) ' 



(4.13) 



where 



(4.14) 



pia;p{A) 



and 



/ gia;p+i(A) 



T^n-p(A) 



I 



(4.15) 
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are p x p and {n — p) x (n — p) matrices respectively and 

1 

.l-A"' ' (1-A) 
Now define a particular matrix solution Q of the Lax pair as 



^iW = a^ ( + ) . (4.16) 



e = (^(/i)Li , ^(/x)L2) , (4.17) 

where Li is an n x p constant matrix of p column vectors and L2 is the orthogonal complementary 
n X (n—p) matrix of (n—p) column vectors. The columns of Li span a p-dimensional subspace U of 
C", and those of L2 span the orthogonal subspace V. The projector P is completely characterized 
by the two subspaces U = ImP and V = KerP given by the condition P-^U = and PV = 0. Let 

us write Li = I ^ 1 and L2 = ( ^ 1 , where A, B, C and D are constant p x p, {n — p) x n, 

px (n — p) and (n — p) x (n — p) constant matrices respectively. Given this, the nxn matric is 
given by 

e=( "-'""-^ V (4.18) 

We now define the projector P in terms of the matrix <I> = ^'(^)Li = (0i, • • • 0p) given by 

= (^i,---,^p) 

\ Wn-MB I ■ 

The projector is thus given by 

\ Wn-p{p)B/^A^w},{ji) Wn-p{^i)B/\B^wl_p{^i) I' 

where A"^ = A^wl{fx)Wp{ii)A + B^wl_p{fi)Wn~p{lj)A. The Darboux matrix D{\) can now be 
constructed to give explicit soliton solution of the SU (n) HM model. To elaborate the result more 
explicitly, we proceed with the example of SU{2) HM model. 

For the SU (2) model, the equations (|4.10p and (|4.11|) become 
Then U"^ becomes 

C/2 = /. (4.21) 
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The Lax pair (jl.lOp - (jl.lip for the SU{2) model can be written as 

d,^ix,t;X) = ^-^Uix,t)^{x,t;\), (4.22) 

dt^ix,t;\) = (^.^-A^U + -^^^UU,^^{x,t;X). (4.23) 

If we take trivial solution (as seed solution), single soliton and multi-soliton solutions can be ob- 
tained by Darboux transformation as explained above. 

We take the seed solution to be 

^0 = ^ = ( -i ) ■ ^^-^^^ 
The corresponding solution of the linear system ()4.22p - ()4.23p can be written as 

Q g V(l-^) (l-A)^ J 

Take Ai = and A2 = fi, the constant matrix A is given by 
and corresponding 2x2 matrix solution B becomes 



_g V(1-m) (1-m)^ / e V(1-m) (1-m)^ / 



The matrix M is given by 

M = GAe"\ 

1 / /ie" + /2e-" (/2-^)e™ \ 
+ e-« (A - /^) e-™ /le" + //e"" J ' 

where the functions u{x,t) and v{x,t) are defined by 

/ ^ /I 1 \ . 1 1 I 
u{x, t) = 1 7:; r — t: r X + 4i ^ ^ t, 

^ ^ V(i-m) (i-/i)y v(i-/i)' {i-fif ' 



(4.28) 
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Let us take the eigenvalue to be // = e^^. The expression ()4.28p then becomes 



^ j cos + i sin tanh u — i (sin 6'sechii) e^'' i (4 30) 

-i (sin Osechu) e~^'" cos 9 — i sin 6 tanh u 



and the corresponding Darboux matrix D (A) in this case is 



I A — cos 9 — i sin 9 tanh u i (sin 9sediu) e^"" \ , , 

D{X)= { _■ ^ ' . (4.31) 

i (sin 9sechu) e A — cos 9 + i sin 9 tanh u 



Comparing the above equation with (j3.12p . we find the fohowing expression for the projector 

p ( 2e"sech« — 2e^''sechti \ 
\ — 2e~^*'sechti 2e~"sechii / 



Using ([32]) and (g^D, we get 

U[l] = { ' 
where 



mn = "t , (4.33) 



C/3 = 1 — (1 + cos 0)sech^ii, 

U+ = U- = -ie^" [(1 + cos 6l)tanhti + i sin 9] sechn. (4.34) 

From equation (j4.34p . we see that C/^[l] = — t^[l] and Tr(C/[l]) = 0. Therefore equation (j4.34p is 
an exphcit expression of the single-sohton solution of the HM model based on SU{2) obtained by 
using Darboux transformation. Similarly one can calculate explicit expression for the multi-soliton 
solution of the model. The expression (j4.34p is similar to the expression of the single soliton given 
in [2]. 



5 Concluding remarks 



In this paper, we have studied GHM model based on general linear Lie group GL{n) and expressed 
the multi-soliton solutions in terms of the quasideterminant using the Darboux transformation 
defined on the solution of the Lax pair. We have also established equivalence between the Darboux 
matrix approach and the Zakharov-Mikhailov's dressing method. In last section we have reduced 
the GHM model into the HM model based on SU{n) and calculated an explicit expression for 
the single-soliton solution. It would be interesting to study the GHM models based on Hermitian 
symmetric spaces. We shall address this problem in a separate work. 
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